Preprint typeset in JHEP style. - HYPER. VERSION hep-th/0104017 

YITP-SB-01-13 

Noncommutative Multisolitons: 

Moduli Spaces, Quantization, Finite 6 Effects and Stability 

Leszek Hadasz 1,2 , Ulf Lindstrom 3 , Martin Rocek 1 , and Rikard von Unge 4 

' C.N. Yang Institute for Theoretical Physics, State University of New York 

O ■ Stony Brook, NY 11194-3840, USA 

2 M. Smoluchowski Institute of Physics, Jagiellonian University 
Reymonta 4, 30-059 Cracow, Poland 

m 



3 Institute of Theoretical Physics, University of Stockholm 
Box 6130 



> 

t^- ' S-113 85 Stockholm, Sweden 

O - 

■^J- ■ ^Institute for Theoretical Physics and Astrophysics 

Faculty of Science, Masaryk University 
Kotldfskd 2, CZ-611 31, Brno, Czech Republic 

leszekOinsti . physics . sunysb . edu 

Qj' 

' ul@physto.se 

rocekOinst i . physics . sunysb . edu 
unge@physics .muni . cz 

% 

Abstract: We find the A^-soliton solution at infinite 6, as well as the metric on 
the moduli space corresponding to spatial displacements of the solitons. We use a 
perturbative expansion to incorporate the leading 0~ x corrections, and find an effective 
short range attraction between solitons. We study the stability of various solutions. 
We discuss the finite 9 corrections to scattering, and find metastable orbits. Upon 
quantization of the two-soliton moduli space, for any finite 6, we find an s-wave bound 
state. 
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1. Introduction 

Recently [Q it was realized that one can construct stable soliton solutions in noncommu- 
tative scalar field theory even though such solitons do not exist in commutative scalar 
theories in higher than two dimensions. The solutions are particularly simple when the 
noncommutativity parameter 9 — > oo, where one finds an infinite dimensional mod- 
uli space of solitons. This program has also been extended to noncommutative gauge 
theories % |, g §, §. 

These solitons have found an application in the context of tachyon condensation 
where D-branes can be found as soliton solutions on higher dimensional non-BPS D- 
branes. By turning on a B-field one makes these non-BPS D-branes noncommutative 
and the soliton configurations studied represent various types of lower dimensional 
D-branes g. 



1 



There also seems to be a place for application of these solitons in a noncommutative 
description of the Quantum Hall Effect |, [H| [11], [12], [R| H§. 



Motivated by these developements, we have studied what happens when one scat- 



ters noncommutative solitons. In [16] we analyzed this questions using moduli space 
techniques, and found a Kahler metric on the moduli space somewhat analogous to the 
metric on the moduli space of two magnetic monopoles. A natural generalization of 



the results in [16] is to find the moduli space metric for iV solitons. In this paper we 
find a simple and elegant expression for the Kahler potential for the general case. 

The analysis in JT| was mainly done at infinite 9 but a program to find corrections 
to the solitons at finite 9 was initiated. This was followed by studies at finite 9, 
both numerically [|TT| , |18|[ and theoretically [[19], 2UJ. This topic is important, since one 



would like to know if the solitons are stable at finite 9, and if they are, how many of 
the infinite number of moduli survive. We study this issue and find that at finite 9, 
nonradial excitations, which were ignored in [|I7|, [18], [19|, ^Oj, destabilize all "excited" 



soliton states, and leave only the basic iV-multisoliton solutions. 

Some quantum issues have also been studied in [F2T], EK|. However, the discussion 



in [[Hj] involves averaging over nonradial modes, which we find play an essential role. 

This paper is organized as follows: In section 2, we construct the general iV-soliton 
solution at infinite 9. We find the metric on the moduli space corresponding to spatial 
displacements of the solitons, and discuss the three-soliton case in detail. In section 
3, we introduce a perturbative expansion that allows us to incorporate the leading 
9~ x corrections. In the two-soliton case we find an effective short range attraction 
between solitons. In section 4 we use these perturbative results to study the stability 
of various solutions. In section 5 we focus on the two-soliton case and discuss the 
finite 9 corrections to scattering. We find a range of interesting phenomena, including 
metastable orbits. In section 6, we quantize on the two-soliton moduli space, and find 
an s-wave bound state for any finite 9. 

While writing of this paper, we became aware that results which have some overlap 
with our results were presented in p3]. 



2. Multisolitons at infinite 9 

2.1 Multisoliton solutions 

The two-soliton solution at infinite 9, constructed in |L], is 

<l> 2 = \(\z + )(z + \ + \z„)(z-\) , (2.1) 
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where A is an extremum of the potential V(4>), 

M = } Z)M ~ Z) , (2.2) 
^(lie-^l 2 ) 

and \±z) = e -2l z l 2±2at |0). This can be generalized to the n soliton case as follows: Let 
z a ,a = 1, ... ,n be pairwise different complex numbers satisfying the center of mass 

n 

condition z a = 0. With 



-\\z a \ 2 +a) z 



|0> , (2.3) 



the multi-soliton solution is 



i = ^E \z a )A${z p \ = \P n , (2.4) 

a,f3=l 

where A is the n x n matrix, 

A afi = (z a | z p ) = e -5d*«-^l 2 +^-^) , (2.5) 

and P n is a rank n projection operator onto the linear subspace of the harmonic oscil- 
lator Hilbert space 7i spanned by the vectors \z a ). 

For large separations, i.e., \z a — z@\ 3> 1 for all a ^ f3, 



|2 



$n= E A l^aX^al+O e-l^l , (2.6) 
a=l ^ ' 

i.e., in this limit, $ n describes n well separated level solitons. To study the limit 
z a — > it is convenient to introduce a new basis: 

| Ul ) = \e Zl ) = |0)+O(e) , 



^- (e^l 2 |ez 2 > - |«i>) = \l) + 0(e) , 



«2 = - f 

l«3> = 17 ^ V ( e|e2 ' 23 ' 2 1^) - |«i> - - \u 2 )) = |2) + 0(e) , 



where e is a small parameter and \i) = (a)^j |0) . P n can now be written as 

n n—1 

P n = Yl M ^/S)" 1 ( U p\ = Yl I*) (*l + C> ( e ) ' 

o,/3=l j=0 
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and hence, in the limit e — > 0, $ n describes n solitons from the up to the n — 1 
harmonic oscillator level, all at the origin. 

In the generic case of different z a , P n is unitarily equivalent to the projector onto 
the subspace spanned by the vectors \i) with i < n. To construct the unitary trans- 
formation explicitly, we diagonalize the matrix A (which is hermitean and positive 
semi definite) . Let uL) and denote its (orthonormalized) eigenvectors and corre- 
sponding eigenvalues: 

^4a/3^/3(7) = 0-( 7 )f Q ( 7 ) , V( a ) ' = 8 a/3 ■ (2-7) 

Defining 



w aP = (2. 



and the Hilbert space vectors 



KHEMw^, (2.9) 
p=i 

we find 

n 

(w a \wp) = 8 a p , and $ n = A ^ \w a ) (w a \ . 

a=l 

For n = 2 this procedure - as expected - gives {|u>i) , |w 2 )} = {|^+) , \ z -)}- 

Finally, we take any orthonormal basis in 7i whose first n vectors coincide with 
\w a ) and denote it by {\wj)}. For 

oo 

tf = £k>0'-i| (2.io) 

we have 

uW = u ] u = i , 

and 

$n = xu (jr \i) ul (2.ii) 

In a completely analogous manner, one may construct excited multisoliton solutions 
from states \z a ,n) = (a') n \z a ) (the excited two-soliton case was worked out in detail 



in 



lqj). However, as shown in section 4, such solitons are all unstable for any nonzero 



9, and hence we do not discuss them further. 
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2.2 Moduli spaces 

The metric on the multi-soliton moduli space is Kahler for any n. Up to a constant 



^ rFml • 

e 








9z a zp 






1 


9z a Zfs — 


A2 




1 


9z a ZfS 


A2 



(2.12) 



Straightforward calculation gives 



9z a Zf3 — 9z a z fj — , 



while 



9z a zfs = A a p A Pa + zpAfaZa - AfrZ^JzsAga = d Za d 2p K (z, 

with 

n 

K (z, z) = ^ \ z a\ 2 + logdet A . 
2.3 The three-soliton case 



(2.13) 



(2.14) 



a=l 



As a simple example, we consider the three-soliton metric in detail; the matrix A 
becomes 

/ l e ~\\zi\ 2 -\\z 2 \ 2 +z lZ2 e -|ki| 2 -5^3| 2 +5i^3 \ 
C.C. 1 

C.C. C.C. 1 



e -^\z2\ 2 -l\z3\ 2 +Z 2 Z 3 



V 



(2.15) 



-zi - z 2 



and gives the Kahler potential (using that z% 

J£ _ J n ^ e 2\zi\ 2 +2\z 2 \ 2 +z 1 Z2+z 2 z 1 _|_ e ~\z 1 \ 2 -\z 2 \ 2 +z 1 z 2 -2z 2 z 1 _|_ e -\zx\ 2 -\z 2 \ 2 -2z~iz 2 +z 2 Z! 

_ -2\zi\ 2 +\z 2 \ 2 -ziz 2 -z 2 z 1 _ \z 1 \ 2 -2\z 2 \ 2 -z 1 z 2 -z 2 zi _ \zi\ 2 + \z 2 \ 2 +2z 1 z 2 +2z 2 zi 



(2.16) 



This has the two-soliton metric as a subspace: If we take one of the solitons (say z%) 
far away and fix its position we can choose the coordinates to be 



zs = ---C, 



(2.17) 
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where ( is the relative distance between the soliton at Z\ and the soliton at z% and is 
taken to be much smaller than z-i- Inserting this in ( [2.16| ) we get 



^|z 2 | 2 + ln(2sinh(2CC)) • (2.18) 

Thus the geometry factorizes into two pieces, one, coordinatized by Z2, which is flat, 
and another, coordinatized by £, with precisely the Kahler potential of the 2 soliton 
moduli space, including the conical singularity when the solitons coincide. 

We can also study what happens when all three solitons come together. We study 
the most symmetrical case where the solitons are at the same distance from the origin 
and separated by the angle ^r. In that case we can choose the coordinates 

Z\ = C, 

z 2 = uj(, (2.19) 

Z3 = ^ 2 C, 

2-n-i 

to = e 3 . 



For small values of ( we get the Kahler potential 

^ = ln(|C| 6 + ^ + ..-) , (2.20) 

giving rise to a conical singularity of the type 

ds 2 = r 4 (dr 2 + r 2 d6 2 ) . (2.21) 
This lead to a scattering angle of 



3. Finite 9: perturbation Theory 

We now consider the finite 9 corrections to the soliton solutions. 

At finite 9 one has to include the derivative terms in the energy functional: 

E = ^Tr (1 [a , 0] [0 , at ] + 9V (0)) . (3.1) 

If we make the ansatz 

= A (p + ia) (3.2) 
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and use the conditions V'(0) = V'(X) = 0, then (|3.1|) gives 1 



E = — Tr 

T 



^V(A)P + y[ a ,P][P,at] + ^B[[P, a t],a] 



(3.3) 



1 1 
+ -V"(X)PBPBP + -y"(0)(l - P)B(1 - P)B(1 - P) 



+ o e 



]-2 



Extremizing with respect to the perturbation B, to leading order we find 



V"{X)PBP + V"{0) (l-P)B(l-P) = [a,[P,a^]} = A P . 



(3.4) 



This implies that P has to satisfy the consistency condition 



[P,A P ] = 0. 



(3.5) 



Fortunately, this relation is fulfilled for all the soliton solutions we consider. 
We can write the solution of (R.4I) in the form 



B 



1 



V"{\) 



PApP 



1 



V"(0) 



;i - P)A P (1 - P) + PX(l - P) + (1 - P)X ] P , (3.6) 



where X is an arbitrary operator; it drops out of both ( |3.4| ) and - as a consequnece of 
(|3.5|) - the energy ( p.3|) to order 6~ l , and hence we choose X = in what follows. 
A special class of projectors that satisfies (|3~5| ) is given by solutions of the equation 



[1 - P)aP = . 



(3.7) 



It includes the projectors that give the n soliton solutions ( |2.4j) . 

It is trivial to see that ( |3.7|) implies (|3.5[) , while the reverse implication does not 
hold - the projector \n){n\ with n > 0, satisfies ( p.5|) but not (|3.7|) . 

For operators P satisfying (|3.7|) , the energy (|3.3j ) can be rewritten in a simpler form 



E = ^Tr ^V(A)P + y p" 



(3. 



vrA 2 



1 1 

+ 



6g 2 \V"(X) V"(0) 



Tr 



2(Pa t a) 2 -2P(a t )V + P 



t ^2^2 



1 No order 2 terms in contribute to the energy at order 9 1 because they are multiplied by 
V'(X) or V (0), which both vanish. 
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3.1 Explicit solutions 

For one-soliton states, our perturbative expansion reproduces the results in Appendix 
A of to the appropriate order in 6~ l . 

For P = \z+) (z+\ + \z— ) (z—\ (|2.1|), the simplifying condition ( |3.7| ) is satisfied. 
One may calculate the first order correction B to $ from 



B = —^—PApP + —^— (l-P)Ap(l-P) , (3.9) 

V"(x) v"(oy J y ' ' y ' 

by substituting A P = [a ,[P,at]\ as above, but to find the corrections to the energy, 
it is easier to use the expression ( |3.9|) , simplify using the cyclicity of the trace, and the 
condition ( [3.7D , and only then substitute the explicit form of P. 
At order 9 the energy is a constant: 

E = -^26V(\). (3.10) 
F 

One might expect that at the lowest non-trivial order the energy could depend on the 
relative position of the solitons z, but this dependence cancels: 

Eh. = %\ 2 . (3.11) 
F 



(For the unstable excited states of [IIJ, E\ = (4n + 2).) The z-independence at the 
lowest nontrivial order implies there is a range of energies for which the moduli space 
is an accurate description even at finite 9. At next order w find 

= vrA 2 ^(A)+^(0) f / 2zz V\ 
2 g 2 6 V"(X)V"(0) \ ^sinh(2^) J J ' 1 ' ' 

The z-dependent part of this produces an attractive force between the solitons. How- 
ever, it is very short-range, vanishing as e~ 4 ' 2 ' 2 . For small \z\, the potential between the 
solitons goes smoothly to a finite constant. (One may consider solitons such that the 
false vacuum is a maximum rather than a minimum: V"(A) < 0. Then the potential 
between the solitons could be repulsive or vanish; however, in this case the solitons are 
unstable.) 

3.2 Corrections to the moduli space metric 

Using the perturbative scheme developed above, we can compute the leading corrections 
to the level zero n soliton metric ( |2.13| ). To order 9' 1 we have: 

9cfl = T^Tr (d a *5 p *) = g% + - Q g% + 0{6- 2 ) , (3.13) 
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where 



9% = \ ( Tr (d a PndpB n ) + Tr (d a B n 8 p P n 



(3.14) 



and d a = d z <* and dp = d 2 /3. For 



(3.15) 



we find: 



0, 



l'"iA) 



taA^zp — A a/3 — ^ A a lz\A\ K z K A K p 

A,K=1 

n 

ZpAp a Z a + A/3 a — A^ZjA^ Z\A\ a 

7,A=1 



(3.16) 



(3.17) 



_ y„,Q\A a p I A^ + 2z f3 A f3a z a + ^ A^zaAa^A^^A^zaA 

v / \ A,k,7,A=1 

n 

- ( 2 + + - z aZp) ApxZ\Allz K A h 

To find the metric on the two-soliton relative moduli space, we set z 1 = z, z 2 
and hence 



Xa 



9zz oc (gn - g X 2) . 



(3.18) 



Explicitly, 



A 



1 e~ 2 W 5 
-2M 2 1 



(3.19) 



and ( |3.17| ) gives 



ijzz 



cothr 2 



sinh r 2 



coth r 2 
^"(0) 



4r z 



sinh 2r 2 \ sinh r 2 



(3.20) 



where r 2 = 2zz. Note that for V"" (A) = V" (0), this correction to / takes exactly the 
same functional form as the original /. 
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4. Stability analysis 



In IJ it was shown that there exists a path in field space interpolating between field 
configurations corresponding to the operators \n)(n\ and |0)(0| and along which the 
gradient energy decreases monotonically. This path is given by |a)(a|, < a < |-, 
where 



\a) = cos a \n) + sin a; |0) 



(4.1) 



Hence for finite 9 the state \n)(n\ decays to the state |0)(0|. 

However, for finite 9 the state |n)(n| differs from the stationary point of the full 
static energy functional by terms of order and the energy of the true solution 
is smaller then the energy corresponding to the |n}(n| state by terms of the same 
order. Due to this energy difference the true level n-soliton cannot decay along the 
|ct)(a| path. Its instability can be, however, still demonstrated with the help of the 
perturbative methods discussed in the section |3|. 

We define 



P a =|a)(a|, A a = [a , [P a , a* 
P a satisfies the consistency condition 

[P a ,A a ]=0 



(4.2) 



(4.3) 



only for a an integer multiple of |, but we can still consider a path in field space given 
by: 



$(a) = A 



P n + \ 



9 \V"(X) 



PAP 

1 a'—^ct 1 a 



1 



7"(0) 



[l-P a )A a (l-P a ) 



(4.4) 



$(0) and $(7r/2) are (up to the terms of order 9 2 ) the true level n- and level 0-solitons, 
and the energy of the $(a) configuration is given by (for n > 1) 



E(a) 



& V(A) + ^ 



g- 



1 + In cos a 
1 



9V"{\) 



2ncos a 



ev"(o) 



1 + 2n cos a + 2n cos a 



+ O (9- 2 ) 



When 9V" 3> n, the function E(a) decreases monotonically for < a < |, and hence 
the \n) (n\ soliton is unstable for n > 1. For n = 1 the resulting formula for E(a) is 
slightly different, but the conclusion is the same. 
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The analysis above can be extended to the case of the solutions corresponding at 
9 = oo to the projection operators 



I C IV . 



(4.5) 



If there is a "gap" in the set I, i.e., for some m < n we have m ^ I and n G I, then 
using the projectors 

-P m ,n(«) = (cos a |n) + sin a |m)) ((n| cos a + (m\ sin a) , (4.6) 

we can construct a "decay path" as in ( |4.4[ ). This shows that potentially stable, radially 
symmetric, level n-solitons <f>r n ) must approach 



too 
\n) 



\Y,\k) (k\ 

k=0 



(4.7) 



as 9 — > oo. We now consider the stability of such states. 
For ijj of the (general) form 

oo 

^({f k },u) = \J2fk^\k){k\u 

k=0 



where U is some unitary operator, the energy functional (|3.1| ) can be written as 

1 



E 



a 2 ^ 



ev(\f k ) + -x 2 f k 



9 \ 2 oo oo 



2ttA 



2 oo oo 



k=0 1=0 



p=0 



with U pJk = (p\U\k). 

The radially symmetric states <f>( n ) have the form 

oo 

0(n) = Cfc \k) (k\ 

k=0 

where, for large enough 9, the perturbative analysis gives 

'1 + (#- 2 ) k < n , 



2(n+l) 

ey"(A) 



+ o(r 



Cfc = < 



0V"(O) " r ^ 



3 — 2>| 



k = n , 
k = n + 1 , 
> n + 1 . 



(4.8) 



(4.9) 



(4.10) 
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To check stability of this solution we consider 



, ) + 5<f ) = XJ2 (c k + e5c k )U^\k) (k\U 

k=0 



(4.11) 



where 5c k are arbitrary real parameters and U = e ieT with arbitrary hermitean T. Using 
we find 



f-J5[0 (B) + S4>] = 9 £) 7(cfc) + A 2 £ (Yfc + -) 4 - (fc + 1) c fc c fc+1 

Z7T k=0 k=0 v v 2/ 

1 oo 

+ 9^ E f^"(Ac fc ) (5c fc ) 2 + (fc + 1) (5c k+1 - 5c k f) 

oo oo r _ 

+ e 2 X 2 E E H / ~ fc ) c fe l T ^| 2 ~ °k c i ^ T Ki-i ~ Vk+lT k+hl 

i — n /— n ^ 



(4.12) 



fc=0 ;=o 



(k + 1) (\T k / + \T k+ltl \ 2 



(k + !)(/ + 1) (T k> iT* 1 + c .c. 



The second line of this equation shows that as long as V"'(Acfc) > for all k, 4>{n) is 
stable against "radial" perturbations c k — > c k + e5c k . 

To check stability against unitary rotations we use the perturbative form of (/>( n ) 
(|4.9|J4.10Q . We denote by £ the T-dependent part of ( |4.12|) . Terms of order 9° in £ can 
be written in the manifestly positive semi-definite form: 



= E E |vTT fc>1 - v / A : r fc _i ii _i 

fc=0 «=n+2 



(4.13) 



In general, these dominate any terms that are lower order in 9; however, (|4.13|) has 
zero-modes, so we need to consider terms of order 9~ 1 : 



£\ 



2(n + l) 
9V"{\) 

2(n+l 
+ ' 



k=0 



I 2 

\ Vn + lT k<n+ i - VkT k -l,n ~ ~ n )\Ti, n \ 

1=0 

oo 2 00 

Y Wn + lT U) i - VlT n+1: i +1 - ~n- l)|T ijTl 



l=n+l 
oo 



9V"(0) 

2(71+1) 

^"(A) /= t 

2(n + l) A I ^— , T rT^ 
g^//( ) E | V n + 2T fc , n+ 2 - VkT k _ hn+l 



(4.14) 



Some of these are positive as they stand, and some are obviously dominated by the 
O(9 ) terms, e.g., the negative sum proportional to yjrjyj together with the terms in 



12 



( 4.13|) with k = n can be written in the form 
which is positive for fy^j < 1- Using the identity 



Vn + lT nj „ + i — \fnT n _i^ n + |T n _ ljn | — |T n)fl+1 | — y/nT n>n +i — \/n + lT n _ 1>n 
we can rewrite the remaining, proportional to V }, X \ terms in (|4.14 ) as 



2(w + l) 
BV"{\) 



V»{\) 
n-l 



y/nT n!n+1 - \/n + lTn-l,n + Y \^ n + lT k,n+l ~ VkT k -l,n 

k=l 

n—2 oo 

+ (n + l)|T , n+1 | 2 + X> ~ fc)|T fe , n | 2 - £ (Z - n)|T, 



-ni 



fc=0 i=n+2 

These terms are potentially dangerous only along the zero mode direction of £\eo with 
some nonvanishing T n j, I > n + 2. Such modes obey the condition 

2 



n oo 



Y \VlT k: i - VkTk-ij-! = . 

fc=0/=n+2 

which gives \Tqi\ = for I > n + 2 and the recursion relation 



a.-,/ 



u-i 



(4.16) 



(4.17) 



this implies 



For Z = 2n + 1 ( P~TTD gives 



T 



for l>2n + 1. 



(n + 1) |T„ j2ri+1 | 2 = 
while for Z = n + p with 2 < p < n we get 



((»+!)!)' 
(2n+ 1)! 



|To. 



n+l 



.77 



(n — p) ! (n + p) ! 



n— p,n | 



and the last three terms in ( |4.15| ) give 



ra-2 



{n + l)\T 0in+1 \ 2 + J2(n-k)\T k , n \ 2 - £ (l-n)\T n>l \ 



k=0 



l=n+2 



[n ■ 



l)(|T , n+1 | 2 -|T n 

,2n+l | 2 ) + YP (\ T n-p,n\ 2 ~ \T n 



p=2 



n+l 



((^ + 1)!) 2 
(2n + 1)! 



\T , n +l\ 2 + 1 



ni. 



p=2 



(n — p) ! (n + p) ! 



|-^n— p,n| ^* • 
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Similar analysis applies to the terms proportional to v ,}^ ■ 

The remaining zero modes of Q4.12|) (apart from the exact, translational zero mode 
which - for 9 —>■ oo - is just T n>n+ i) are given by T^i with k > n + 1 or with I < n 
(thanks to the hermiticity of T we can always choose k to be smaller than I). The first 
set of these zero modes (with k > n) is irrelevant - they correspond to the unitary 
transformations that do not change the state to order 9~ 1 . 

The other set of zero modes (with k < I < n) corresponds to rotations that act on 
the nonzero c^'s and do change the soliton. To check if they can destabilize the soliton 
we would have to extend the perturbative analysis to higher orders in 9~ l . Fortunately, 
they are absent for the most interesting case, that of two solitons (i.e. for n = 1). 

5. The geodesic equation at finite 9 

We now study classical scattering in the presence of a generated potential U. 

5.1 Integrating the classical equations of motion 

For concreteness and simplicity we consider the two-soliton case (|3.12|) : 

U(r) = - (±] V " {X) + r/(0) 2 (5 1) 

K> \29J V"{X)V"{0) \smh{r 2 )J ' y ' } 

The effective action for the two-soliton system is 

2vrA 2 r , (9 



9 2 J \2 



dt[~{g r y + g M r)-U{r)\ , (5.2) 



where the metric g is given by ds = f(r)(dr + r dd 1 ) and / is as in [16 



/(r) = coth(r 2 ) - , ^ . (5.3) 
smh [r z ) 



In principle, we should consider the corrections (|3.20|) to the moduli space metric; as 



discussed below, for large 9, these can be ignored. In contrast, the potential U is 
important even though it is also suppressed at large 9. 
Varying the action leads to the equations of motion 



^r: r r 2 + r;/ + ^ = o, (5.4) 

+ 2T%H) = . (5.5) 
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The second equation is the same as for the case without U and the solution is 

= ~2~cT~\ ( 5 - 6 ) 
r z j[r) 

where / is an integration constant corresponding to the angular momentum. The first 
equation can then be written as 

An integrating factor for this equation is rf from which one finds the solution 

= 2^ - 41 , (5.8) 
where E is an integration constant with the interpretation of the total energy of the 



system. In the same way as in the case without a potential [fLq| , this leads to scattering 
trajectories found from the integral 

/r (Jo 
1 2 (5-9) 

here b 2 = is the impact parameter as in |16| (the 9 dependence arises because we 
have rescaled the coordinates). The finite 9 correction to the geodesic scattering picture 
can therefore be found by using a corrected function 

/» = f(r) (l - ^) • (5.10) 

Since U is attractive (negative) there are no extra divergencies in the effective / as 
compared to /. If U had been repulsive, but of the same functional form 2 , it would 
have made the effective / more repulsive. 

We can make some estimates of the validity of our approximations by restoring the 
dimensions of the coordinates: r — > Since E = ^{^v 2 + U), which corresponds 
to a particle of effective mass moving in a potential U, we can find a range of 
velocities where the moduli space approximation should be good. For the correction to 
the classical result in ( |5~9l) to be small we need ^ <C 1 leading to 

1 v"W + v"(Q) , n) 

V ^ 29 V»(\)V»(0) ■ {b } 



2 This could happen for a potential V where the false vacuum corresponds to a maximum; however, 
such solitons are unstable. 
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However, for the adiabatic approximation to be valid, the momentum transfer must 
remain sufficiently small so that fluctuations out of the moduli space are suppressed. 
In our case there are several possiblities since potentials for different fluctuations appear 
at different orders in perturbation theory. Even if we do not have the exact potentials, 
we can estimate their strength from the general behavior of perturbation theory. Most 
fluctuations have potentials already at 9 = oo. They are not excited as long as 

^ « ^ (5.12) 

Other fluctuations get a potential only at first order in perturbation theory. They are 
not excited as long as 

v 2 < 1 ; (5.13) 

this simply means that the motion remains nonrelativistic. For the two-soliton case, 
we have checked that the fluctuations with lower energies correspond to motions of 
the solitons, which we do not want to restrict. Higher soliton scattering requires a 
higher-order analysis. 

5.2 Trajectories 

It is interesting to investigate some explicit cases for the scattering trajectories of the 
previous section. We have prepared movie clips in MPEG format 3 . The | first movie 



fchows the behavior for large values of the impact parameter. The solitons just pass each 



other with no scattering taking place. In the second movie" the right angle scattering 
for small impact parameter b is shown 4 . Notice that this qualitative behavior is true 
irrespective of the value of the total energy E since it only depends on the value of 
the function / at large or small r. However, in the presence of the attractive potential 
U(r) and for small enough energy > 3.86) we find new qualitative behavior shown 



in the | third movie .] We get a metastable orbit where the solitons circle around each 
other for some time before they scatter to infinity. These results are summarized in the 
following picture where the exit angle is plotted as a function of the impact parameter 
in the case where ^ = 5. 



3 If the reader's viewer does not support hypertex, the three movies can be found 
at jittp : //www . physto . se/~unge/traj 1 .mpg , http://www.physto.se/~unge/traj2.mpg, and 



tittp: //www. physto . se/~unge/traj3 .mpg 



4 The second movie has time slowed down by a factor of 1000; to keep the file a manageable size, a 
smaller spatial region is shown. 
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Figure 1: The exit angle in degrees as a function of the impact parameter for ^ 



One may also ask what happens when we include the corrections to the moduli 
space metric ( P^D|) . When V"{\) = V"{0), f is rescaled by (1 + f ), which can be 
absorbed by a redefinition of b. However, when V"(X) and V"(0) are not equal we find 
two different behaviors: If V"(A) > V"(0), the trapped orbit effect above is suppressed, 
that is, it appears for smaller E, or larger than before. On the other hand, if 



V"(0) > V"(A) we get an enhancement of the effect, 
even for Uq = if V"(A) is small enough! 



In fact, trapped orbits appear 



6. Quantization 

In this section we discuss the quantization of the effective hamiltonian that describes 
the motion of solitons. We focus on the two-soliton case. The Schrodinger equation for 
this problem can be written as 5 

(-V 2 + U) tjj = Etfj , (6.1) 

where 

V 2 = ^d r (rd r ) + -±- 2 dl , (6.2) 
J r J r 



5 There is really a factor of in front of the V 2 operator but we can soak it up in a redefinition 
of U and E so that they become dimensionless. 
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and f(r) is the metric (|5.3| ). In the absence of a potential, this Hamiltonian operator is 
positive 6 ; consequently, it cannot have any bound states. For any attractive potential, 
no matter how small, an s-wave bound state appears. Thus the potential that we found 
perturbatively (|3.12|) induces such a bound state for any finite value of 9. The potential 



is rotationally symmetric and hence for ip = x(r)e* Zl? the equation reduces to: 



~b(ra r ) + j pi + u)x = Ex. (6.3) 



We may estimate the energy of the bound state as follows: for sufficiently small r, 
/ « |r 2 and U ~ Uq < 0. This admits a solution with energy E of the form of a Bessel 
function J Q (kr 2 ) for k 2 = (Uq — E)/6. For large r, / ~ 1 and U ~ 0; this admits a 
solution of the form A Ko(y—E r), where A is some normalization constant and Kq is 
a Bessel function of the second kind. Matching these at some intermediate value r m 
where U(r m ) — > E gives 



2kr m J 1 (kr 2 m ) y/^EK^y/^E: 



J (kr 2 J K {V=E' 
for small enough \Uq\, \E\, we find 



(6.4) 



E « - J-e-VC^) (6 5) 

Of course, the details of f(r) and t/(r) correct the solution and the energy, but they 
cannot change the qualitative behavior. 

Because of this s-wave bound state, for sufficiently small soliton energy the cross 
section becomes very large; this ruins the moduli space picture. The conclusion is 
therefore the same as in the previous section: for the moduli space picture to be a good 
approximation we need energies in an intermediate range, not big enough to ruin the 
adiabatic approximation but not so small as to see the bound state. 

The question of how the classical metastable states found in the previous section 
appear in the quantum treatment is left to future work. 
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Note added: After completing our work, |24| appeared. This work has significant 
overlap with ours. The authors observe that the solitons naturally obey Bose-Einstein 
statistics, and that imposing this on the solitons replaces the moduli space with its 
quotient by the symmetric group; as pointed out in ||16[ (where, however, we did not 
impose Bose-Einstein statistics), this quotient smoothes the singularities discussed in 
section 2. In the quantum analysis of section 6, Bose-Einstein statistics require us 
to consider only even angular momentum I. Finally, we want to thank the referee 



for making us aware of ||25|| , where noncommutative solitons are considered from a 
somewhat different perspective. 
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